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Abstract 

In this paper, we investigate the Hawking radiation of massive spin-1 particles from 4-dimensional 
Kerr and Kerr-Newman black holes. By applying the Hamilton-Jacobi ansatz and the WKB 
approximation to the field equations of the massive bosons in Kerr and Kerr-Newman space-time, 
the quantum tunneling method is successfully implemented. As a result, we obtain the tunneling 
rate of the emitted vector particles and recover the standard Hawking temperature of both the two 
black holes. 
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I. INTRODUCTION 


Hawking proved that black hole can release radiation thermodynamically, due to quantum 
vacuum fluctuation effects near the event horizon [1, 2], Since then, several early methods 
were proposed to derive the Hawking radiation [3-8], mostly relying on quantum held theory 
on a flxed background. Later, a semiclassical derivation of Hawking radiation as a tunneling 
process, including null geodesic method and Hamilton-Jacobi method, has been developed 
and has attracted a lot of attention. The tunneling methods provided people an alternative 
way of understanding black hole radiation. Both approaches to tunneling used the fact that 
the tunneling probability for the classically forbidden trajectory from inside to outside the 
horizon is given by T = exp (—2Im//h), where I is the classical action of the trajectory. The 
difference between these two methods consists in how the imaginary part of the classical 
action is calculated. 

Null geodesic method was first put forward by Kraus and Wilczek [9, 10] and then de¬ 
veloped by Parikh and Wilczek [11-13]. For this method, the only part of the action that 
contributes an imaginary term is where Pr is the momentum of the emitted null 

s-wave. One can compute the imaginary part of the action by using Hamiltons equation and 
knowledge of the null geodesics. Works based on null geodesic method consist of [14-18]. 

Hamilton-Jacobi method was proposed by Angheben et al . [19], as an extension of Pad- 
manabhan’s works [20, 21]. This method applies the WKB approximation to the Klein- 
Gordon equation, with the Hamilton-Jacobi equation as the lowest order. Then according 
to the symmetry of the metric, one can pick an appropriate ansatz for the action and put 
it into the Hamilton-Jacobi equation to solve. Applying this method, people studied the 
Hawking radiation of various space-times and made a great deal of successes [22-25]. In 
a way similar to Hamilton-Jacobi method, R. Kerner and R.B. Mann [26, 27] applied the 
WKB approximation to the Dirac equation to calculate Dirac particles’ Hawking radiation. 

Kerr and Kerr-Newman metrics describe the geometry of vacuum space-time around a 
rotating uncharged and charged axially-symmetric black holes, both of which are of great 
astrophysical significance or mathematical interest. The Hawking radiations from these 
two kinds of black holes have been studied intensively. Zhang et al . [28, 29] and Jiang 
et al . [30] studied massless particles tunneling from Kerr black hole and that of charged 
particles from the Kerr-Newman black hole in the Parikh-Wilczek tunneling framework. 
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Kernel et al . [27], Li et al . [31] and Jian et al . [32] investigate the tunneling of fermions 
from the Kerr and Kerr-Newman black holes by applying the WKB approximation to the 
Dirac equation. Umetsu [33, 34] studied the tunneling mechanism in Kerr-Newman black 
hole by introducing the technique of the dimensional reduction near the horizon. Chen et 

[35] used the Landauer transport model and one-dimensional quantum channel theory 
to investigate Hawking radiation of photons and fermions in Kerr and Kerr-Newman black 
holes. Vieira et al . [36] and Konoplya et al . [37] precisely studied the Hawking radiation 
of charged massive scalar particles from the Kerr-Newman black hole by analytical and 
numerical methods respectively. All the literature above studied the Hawking radiation 
from Kerr and Kerr-Newman black holes by investigating the tunneling of photons, scalar 
or fermion particles. 

In this paper, we aim to work out the Hawing radiation of massive vector particles from 
the Kerr and Kerr-Newman black holes. Vector particles (spin-1 bosons) such as Z and 
bosons are well-known and play very important role in Standard Model. When we consider 
the massive bosons in Kerr space-time or the uncharged massive bosons in Kerr-Newman 
space-time, the motion of the boson held can be depicted by the Proca equation. In this 
case, we can study the tunneling mechanism by applying the Hamilton-Jacobi ansatz and 
WKB approximation to the Proca equation in Kerr or Kerr-Newman space-time, along the 
same path as in [38-43]. However, the motion of the charged massive bosons in the Kerr- 
Newman black hole should be more complicated than the Proca equation because of the 
nontrivial interaction between the charged massive boson (kP^) helds and the electromag¬ 
netic held. Firstly, we derive the held equation of IK-boson from the Lagrangian given by 
the Glashow-Weinberg-Salam model. Then we apply the Hamilton-Jacobi ansatz and WKB 
approximation to the derived equation in Kerr-Newman space-time. By setting the deter¬ 
minant of the derived coefhcient matrix to zero, a set of linear equations can be solved for 
the radial function. Consequently, we compute the tunneling rate of the charged massive 
vector particles from the Kerr-Newman black hole and recover the corresponding Hawking 
temperature. 

The remainders of this paper are outlined as follows. In Sec. H, taking the rotation ehect 
into account, we investigate the Hawking radiation of vector particles from the Kerr black 
hole and recover the Hawking temperature. Extending this work to the charged rotating 
space-time, both the uncharged and charged bosons tunneling from the Kerr-Newman black 
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hole are studied in Sec. III. Sec. IV contains some discussion and conclusion. 


II. QUANTUM TUNNELING OF MASSIVE VECTOR PARTICLES FROM 
KERR BLACK HOLE 


In this section, we investigate the Hawking radiation of massive spin-1 particles from 
rotating (Kerr) black hole. The line element within Kerr space-time is given by 

2 Mr, 


ds^ = —(1 — )dt^ + ^dr'^ + p^dO'^ 

A 

, 9 9, 2Mra^sm‘^9] , AMrasin^O , , 

(r^ -h a^) H-5- sm^Odp^ -5- dtdp, 






( 1 ) 


where p^ = + a^cos^P, A = — 2Mr + with M as the black hole mass and a as 

the angular momentum per unit mass. To make the event horizon coincide with the inhnite 
red-shift surface, we introduce a new coordinate y = 9 ? —Qt with Q = under 

which the metric ( 1 ) becomes 


ds^ 


Ap 2 


E(r,P) 


dt^ 


^dr'^ -|- p^dO"^ + 
A 


S(r,P) 

p 2 


sin^Phy^ 


( 2 ) 


where S(r, P) = (r^ -|- a^Y — Aa^sin^P. 

In a curved space-time without electromagnetic held, the motion of massive spin-1 vector 
held is depicted by the Proca equation 


Til 

= 0 , 


( 3 ) 


in which 

(4) 

Because of the anti-symmetry of tensor, the Proca equation is equivalent to 

1 _ T77^ 

= 0 . (5) 

By substituting the metric (2) into (5), the Proca equation can be reformulated as a 
tetrad 
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+^3 

<90 
-ds 
do 

-do 

do 

-ft 




P" 


+ ft 


-^—^(ft^3 — ft^o) 


Asin^P 

\/^S 


P 

Ssin^P 

\/^S 


(ft^o — ft^i) 


- ft 


(ft^3 - ft^l) 


|^^(ft^ 2 -ft^ 0 ) 

Ap 2 n 2 ° 

^^^ (ft ^2 -ft^l) 


Ap4 

■x/=^ 


Esin^P 

^/=^ 


(ft^O — ft ^2 
(ft^3 — ft^ 2 ) 


Asin^P 

WEL 

Ssin^P 


(ft^o — ft'J's) 


(ft ^2 — ft^3) 


ft 

+ 

-ft 

+ 


p'^h? 

ft=pA 

y/—gm‘^ 


-^1 = 0 , 


(ft^'l - ft^' 2 ) 


p 2 /j ,2 

'EEmEl 

.Esin^P 

Esin^P/i^ 


^2 = 0 , 

(^3^1 - ft^s) 
^3 = 0 . 


According to the WKB approximation, \1/^ is of the form 


'I'm = C'M(^,A^,x)exp 


■^s{t,r,e,x) 


where S is dehned as 


( 6 ) 

(7) 

( 8 ) 


(9) 


( 10 ) 


S{t, r, 9, x) = So{t, r, 9, x) + hSi{t, r, 9, x) + r, P, x) H-• (H) 

Substituting Eqs. (10), (11) into the Proca tetrad (6)-(9) and keeping only the lowest order 
in h, we get equations of the coefhcients Ca 

^ [C'o(ft^o)' - C'i(ftAo)(ft^o)] + ^ [C'o(ft^o)' - C'2(ft^o)(ft^o)] 

+ [Co(ft^o)' - C3(a^^o)(ft^o)] - Com^ = 0, (12) 

^ [C'i(ft^o)' - C'o(ft^o)(ft^o)] + ^ [C'i(ftAo )2 - C'2(ft^o)(ft^o)] 

[C'ftftAo)' - C'3(a;,Ao)(ft^o)] - Cim^ = 0, (13) 

0 [C'2(ft^o)' - C'o(ft^o)(ft5o)] + ^ [C2{dtSof - C,idrSo){deSo)] 

[C2id^Sof - C3(a^^o)(ft^o)] - = 0, (14) 

^ [C3(ft^o)' - Co(ft^o)(5x^o)] + ^ [C^i(ft^o)' - Ci(ft^o)(9x^o)] 

+ ^ [ 773 ( 5 ,^ 0 )' - C'2(ft^o)(ax‘^o)] - Com^ = 0. (15) 
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Considering the properties of the Kerr space-time, we carry ont separation of variables 


as 


S-o =-Et + W{r)+jip + Q{e) 

= -(E - ]a)t + W(r) + ix + W), 


(16) 


where E and j denote the energy and angnlar momentum of the emitted particle respectively. 

Inserting Eq. (16) into Eqs. (12)-(15), one can obtain a matrix equation K (Co, Ci, C 2 , Cs)^ = 

0 and the elements of K are expressed as 

, J," , /\W'(E-ia) 

^-1-^ - m , Ki2 — -^-, 


Kn = 
K22 = 

K24 = 


Ssin"^^ 


Je{E-jn) 

s(E-jn) 2 ^V^ rp 




P^jiE-jn) -EW'{E-jn) 

y -f^l4 — - ^ ■ 2/1 - 1 -^21 — 


Ap 2 p^ 

zAe:, K 22 =_ 

Esin^P Ap 2 


Esin^P 
,'2„2 

- iCo-i = 


Ap 2 

-JeW 


Esin^P 


-EJeiE-jQ) 


, Ks2 = 


-AJ 01 E' 

9 ’ 

p2 


(17) 


-EiE-jnf , AlE'^ , pY 

2^33 — -;—- 1 --h 


Ap2 


P" 


Esin^P 


2 „ -p^jJe 
m , K 34 = 


Esin 2 p 


2/3 ’ 


-Ej{E-jQ) -AjW -Jej 

-^41 — -X—5-1 -^42 — -5 5 -^43 — ^• 


i^44 = 


Ap 2 

-E(E-jfl )2 


P" 

AIE'^ Je^ 2 

+-^ + ^ - mC 

p^ p^ 


P" 


Ap 2 

where Jq is identified as OoSq. 

The determination of the coefficient matrix should be equal to zero to ensure that Eqs. 
(12)-(15) possess nontrivial solution. By solving detiC = 0, we obtain 

^,2 p^sin^PE2(E — jVtY -\- p^sin^Pm2AE — p2sin2pjg2/^y; _ p^j'^A 


IE'" = 


(18) 


A 2 p 2 sin 2 ^y; 

which means that there are two values contrary to each other that can be identified as the 
radial derivatives of the action Sq. Correspondingly, we have 


IE±(r) = ± 


I p2sin^PE2(E — jT2)2 -|- p"^sin^Pm2AE — p2sin2pjg2/\y; _ p6j2^ 


dr. (19) 


A 2 p 2 sin 2 ^y; 

We mention that 1E+ denotes the radial function of the outgoing particles and 1E_ of the 
ingoing particles. Integrating around the pole at the outer horizon = M + y/— a?., 
we obtain 

H4(r) = ±„(l±!±AAxA±), (20) 


2(r+ — M) 
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where ^2+ is the value of hi corresponding to r = r+. 

So the tunneling probability of the vector particles is 


r = 


Pr 


outgoing 


Pi 
= exp 


exp [—|(Imhh+ + Im0)] 

ingoing OXp [-\{lmW_ + Im0)] 
27rr+^ + 


= exp 


4 

n 


-{E-jn+) 


h r+ — M 

Setting h = 1, then Hawking temperature of the Kerr black hole is recovered as 

1 — M 


( 21 ) 


( 22 ) 


27r r +2 + ’ 

which is exactly the result obtained by other methods. This implies that when black hole 
radiates scalar particle, Dirac particles and vector particles, their tunneling probability and 
Hawking temperature are the same and are not related to the kind of particles. The Hawking 
radiation spectrum of the bosonic particles from Kerr black hole can be deduced following 
the standard arguments [3, 44] 

1 


N{E,j) = 


E-jn^ 

e ^b'^h — 1 


(23) 


III. QUANTUM TUNNELING OF MASSIVE VECTOR PARTICLES FROM 
KERR-NEWMAN BLACK HOLE 


We can obtain the metric of Kerr-Newman space-time by replacing 2Mr in Kerr space- 
time (1) with 2Mr - Q‘^ 

12 ) 


ds^ = -(1 - 


2Mr — , 2 , 2 2 io 2 2{2Mr — Q‘^)asm.^6 

+ ^dr^ + p^dO^ - ^ - dtd^ 


p- 




+ 


2 2 \ {‘2Mr — (5^)a^sin^P 


{P + a^) + 


P" 




with the electromagnetic potential 

= Atdt -\- Aipdp = 


Qr 


, QraAv?6 , 
-dt --vrdP) 


(24) 


(25) 


-|- a^cos^P -|- a^cos^P 

where Q is the charge of the black hole, and A = r‘^—2Mr+a?+Q‘^. Similar to the Kerr black 
hole case, we introduce coordinate transformation y = P — where Q = , , o„ , 

’ ^ -r 1 (r-2+a2)2-Aa2sin^e ’ 

then the metric (24) becomes 

Ap 2 


ds = 


a , P .J„2 , 2 ./)2 , ^)„. 2/1 ._2 


E(r,P) 


dP + ^dP + p"dP" + ^ 
A p 2 


-sin Pdx , 


(26) 
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with the corresponding potential 


— Atdt + A^dx 


Qr{r‘^ + 

S(r, 0 ) 


Qrasin^O ^ 
r2 + a2cos20^^’ 


(27) 


where S(r, 9) = (r^ + — Aasin^6'. 

The charged bosons {W^) behave differently from the uncharged bosons {Z) in the elec¬ 
tromagnetic held. So it’s reasonable to investigate their tunneling processes separately. 


A. Uncharged massive vector particles 


The motion of uncharged massive vector held in the Kerr-Newman background is still 
depicted by the Proca equation (5), except that the metric refers to (26) now. Because of the 
similarity of the Kerr-Newman metric (26) with the Kerr metric (2), it should be trivial to 
extend the study from Kerr case to Kerr-Newman case. One can repeat the process just by 
replacing A, 0 and S with A, 0 and S respectively. As a result, the tunneling probability 
is 


T = exp 


-27r 


ri -I- 


{E-jn+) 


rj_ — M 

and the Hawking temperature of the Kerr-Newman black hole is recovered as 


(28) 


1 r+ — M 
27r -|- a2 ’ 


(29) 


where h= 1 has been assumed, = M+ M‘^ — Q'^ — a? is the outer event horizon of Kerr- 
Newman black hole, and 0+ is the value of O at r = r+. Again the Hawking temperature 
we obtained in (29) is fully in consistence with that obtained by other methods. 


B. Charged massive vector particles 

According to the Glashow-Weinberg-Salam model [45], the Lagrangian of the hP-bosons 
in a background electromagnetic field is of the form 

1 m ^ 

C = -- {D+W+ - D+W+) 

^ ! L 

-^eF>^^W+W-, (30) 

where ± with as the geometrically covariant derivative and A^ as the 

electromagnetic potential of the black hole, e denotes the charge of the hP’*' boson, and 










The equation of motion of the W-boson held can be derived from 

Eq. (30) 

± = 0. (31) 

Since the equation of motion of the boson is similar to that of the W~ boson, they 
should share similar tunneling processes too. For simplicity, we will study in detail only 
the boson case, after which the results can be extended to W~ boson. Because of the 
diagonalization of the metric (26), the equations of motion of held can be reformulated 
as 


1 


v 9 \_9fj.fj.9i 


^ ^ d^W+ - a,lE+ + ^e2l^lE+ - ^e2l,lE+ 


h 




^9^J,^J.9v 


h 


h 


-fE+ + -eF^^W^ = 0, for z/ = 0,1, 2, 3, 


(32) 


where only the summation over p is implemented, and the relation V^lEj^ — = 

d/j^Wif — di^W^ has been used. 

According to the WKB approximation, W^f is of the form 



Cf,{t,r,9,x)ew 


i 

h 


S{t,r,9,x) , 


(33) 


where S is dehned as in Eq. (11). Substituting Eq. (33) and (11) in Eq. (32), the equations 
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to the leading order in h are 


^ [C'o(a,^o)' - C'i(a,5'o)(at^o + eAt)] - 
[C'o(aeS'o)2 - C2{deS,){dtSo + e^)] 

+ ~ ^ 2/1 [^o(^x‘^o + ~ C^{d-^SQ + eA^{dtSQ + eAt)\ = 0, (34) 

Esin 6* 

[C'i(c?fS'o + eAtY — Co{dtSo + ey4i)(c?j.*S'o)] — Cirn^ 

Ap^ 

+^[C,ideSof-C2{deSo)idrSo)] 

[Ciid^So + eA^f - Csid^So + eA^){drSo)] = 0, (35) 

Esin 6* 

^ [C2{dtSo + eAtf - CoidtSo + eAt){dgS,)] - C2mlr 
Ap^ 

[C2{dtSo + eAtf - Ci{drS,){dgS,)\ 

+ ~ ^ 2/1 [^2(^x‘^o + ~ C^{d-^SQ + ey4^)(c?6iS'o)] = 0, (36) 

Esin 6* 

~ [C'3('9fS'o + eAtY — Coi^dtSo + eAt){d^So + eA^)] — 

Ap 2 

+— [C'i(c?tS'o + eAtY — C'i(5r.5'o)(5^S'o + eA^)] 

+— [C'3(c?eS'o)^ — C2{d0SQ){d^SQ + eA^)] = 0. (37) 

Note that the nontrivial interaction term in Eq. (31) doesn’t contribute to the 

motion of the hh-boson held at the leading order in h. 

Considering the properties of the Kerr-Newman space-time, Sq exists with the same form 
as in the Kerr space-time 


5-0 =-Et + W{r)+]ip + Q{e) 

= -{E-3n)t + W{r)+3x + Q{e). 

Inserting Eq. (38) into Eqs. (34)-(37), one can obtain a matrix equation 

K(Co,Ci,C 2 ,C 3 )^ = 0 , 


(38) 


(39) 
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where K is a 4x4 matrix whose elements are 




P 


-2" + 


Ssin ^0 




EiE-eA-jnf , V , (j + eAjV 


/^22 = 
i^33 = 
^^44 = 


H-^ + 


Ap2 Ssin^6* 


— m 


w 


-EjE-eA-jnf Ahh^^ , p\j+eA^) 

Ap 2 p 2 

-E{E - eAt - jQ)^ AW'^ 


EsinA 




Ap 2 


P^ 


ii 

p 2 


m 


Wy 


AW'(E-eA-jQ) _ME-eA-jn) 

-ft 12 — - o -) -^13 — --• 


P" 


P" 


(40) 


p^(j + eA^)(E-eA-jn) -EW'(E - eA - jQ) 

A 14 — -^---, A 21 — -^-. 


Ssin^P 

_ -J,W' _ -f?U + eA^)W' 

-^23 — - 5 -) -^24 — 


Ap 2 


/t'31 = 


A ^ Ssin^P 

-EJq{E - eAt-iA 

—, -ft 32 — 


Ap2 


-AJelh"' 

9 ’ 

p 2 


rx _ -P^(j + eA^)j 0 _ -S(j + eA^){E - eAt - 

A 34 —-~ . n „ -) 2141 —-^-. 


Ssin^P 
-A(j -I- eA^)W' 

-^42 — - o -) -^43 — 


Ap 2 

-Je(j+ eA^) 


Eq. (39) possesses nontrivial solution if the determinant of the matrix K equals zero. By 
solving detJ’C = 0, we obtain 

^,2 p^sin^PE^(E — eAt — jhl)^ + p'^sin^Pm^AE — p^sin^PJ^^AS — p®(j + eA^^A 

A2p2sin^PE 

(41) 

Integrating AW' around the pole at the outer horizon r + = M + obtain 

(42) 


.x. / s , • (^ + «^)(^ - e^f+ - i^+) 
iy±(r) = ±*7r —— 


2(f+ - M) 

where 12 + and At+ are the values of 12 and At corresponding to r = ?+ respectively, iy+ 
denotes the radial function of the outgoing particles and W- of the ingoing particles. Then 
we can obtain the tunneling probability of the bosons 


r(lE+) = exp 


4 

--ImlE+ 

n 


= exp 


27r A + A 
h A — M 


{E - eAt+ - jl2+) 


(43) 
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and Hawking temperature at the outer event horizon of Kerr-Newman black hole is recovered 


as 


T„ = 


1 r+ — M 


^ 27r ’ 

which is fully in consistence with Eq. (29). 


(44) 


The W~ boson case proceeds in a manner fully analogous to the W~^ boson case discussed 
above. Similarly, the tunneling probability of W~ bosons can be obtained as 

27r rl + 


r(lE-) = exp 


:{E + eAt+ — jT2_|_) 


h — M 

and the Hawking temperature in Eq. (44) is recovered once again. 


(45) 


IV. CONCLUSIONS AND DISCUSSIONS 

In this paper, by applying the WKB approximation and the Hamilton-Jacobi ansatz to 
the massive vector held equations, we have investigated the tunneling mechanism of massive 
bosons from the Kerr black hole and Kerr-Newman black hole. The expected Hawking 
temperatures have been recovered and proven to be consistent with black hole universality. 
We should mention that in our investigation, the back reaction of the emitted particle on 
the black hole geometry and self-gravitational interaction are reasonably neglected, and the 
derived Hawking temperatures are only the leading terms. One can compute corrections to 
the massive bosons tunneling probability and the emission temperature by fully taking into 
account conservation of energy and charge. 
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